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Abstract 

We present a model-independent analysis of CP violation, inspired by recent experimental 
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I. INTRODUCTION 


Numerous studies of CP violation have been carried out in quark flavor physics. The b- 
and s-quark sectors have provided a fertile testing ground for the the standard model (SM) 
explanation of CP violation through the Cabibbo-Kobayashi-Maskawa (CKM) matrix. B 
factory experiments played a central role to confirm the CKM framework and to determine 
each matrix element. 

In the charm sector, there are stark differences for the study of CP violation. Mixing 
occurs at extremely small rate in neutral charmed mesons, compared to that in neutral K 
and B mesons. It results in the mixing-induced indirect CP asymmetry being negligible. For 
direct CP violation, the asymmetry can vary greatly depending on particular charmed meson 
decay modes. The largest direct CP asymmetries are expected in singly Cabibbo-suppressed 
(SCS) decays, such as —)■ 7r’''7r“ and —)■ K~^K~, where interference between penguin 

and tree contributions can be substantial. Naturally both theoretical and experimental 
interest has been focused on this type of charm decays. In Cabibbo-favored (CF) decays, 
such as —)■ K~'k^ and —)■ , the favored tree contribution dominates so that 

direct CP violation is negligible. In doubly Cabibbo-suppressed (DCS) decays, such as 

—)■ K^'k~ and , direct CP asymmetries are expected to be negligible in the 

SM, but non-negligible in certain new physics (NP) models. 

Three years ago, LHCb had observed indications of direct CP asymmetry at 3.5cr in 
Aocp K^K~) — acp{D^ —)■ 7r’''7r“) [l|. CDF and Belle also reported similar 

I 0.147)%. These results attracted 


23|, both in the context of the 


results [ 2 I, af, giving a world average of Aacp = (—0.678 
great attention and led to numerous theoretical works 
SM and of various models of NP. For a long time, direct CP violation in charm decays 
was expected to be below 0(10“^). However, some of the recent studies indicate that 
asymmetries at 0(10“^) in these final states may not be excluded within the SM. Since 
then, LHCb has reported new results from which the world average has moved much closer 
to zero: Aocp = (—0.253±0.104)% 24|. The direct CP asymmetries, reported by LHCb 25|, 
for 7r+7r“ and —)■ K^K~ are acp(vr+7r“) = (—0.20 ± 0.19(stat) ± 0.10(syst))% and 

acp{K^K~) = (—0.06 ± 0.15(stat) ± 0.10(syst))%. More precise measurements from future 
experiments, such as upgraded LHCb, BESHI and Belle H, are crucial to provide necessary 
information. Besides CP asymmetries, the very different experimental values of branching 
fractions B{D^ K^K~) and B{D^ —)■ 7r’''7r“) have bee^a long-standing puzzle. Some 


efforts have been made to resolve the issue within the SM 


23 


26|. 
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For hadronic charmed meson decays, there is still no proper theoretical description of the 
underlying mechanism based on QCD. As well known, it is mostly because of the charm quark 
mass of order 1.5 GeV, which is not heavy enough to apply for a heavy quark expansion, and 
not sufficiently light to allow for the application of a chiral expansion. Thus, for hadronic 


D decays, it is not reliable to use the QCD-inspired theoretical approaches satis: 


worked for B meson decays, such as t 
perturbative QCD (pQCD) approach 


le QCD factorization (QCDF) approach 271. 


29 


actorily 
23 !. the 


m and the soft-collinear effective theory j3l[ |. 

ecays of charmed mesons, a 


32 


34| , has been developed. 


Since there is no reliable theoretical framework for hadronic c 
model-independent method, called the quark diagram approach 
In this approach all the decay amplitudes are decomposed into the so-called topological am¬ 
plitudes corresponding to the different topological quark diagrams. Based on flavor SU(3) 
symmetry, the heavy meson decay amplitudes also can be decomposed into linear combi¬ 
nations of the SU(3) amplitudes which are SU(3) reduced matrix elements 35l442l| . This 
approach is equivalent to the quark diagram approach when flavor SU(3) symmetry is im¬ 
posed to the latter. Since each topological quark diagram includes all possible strong inter¬ 
actions to all orders, analyses of topological diagrams can provide information on final-state 
interactions (FSIs). In this model-independent analysis, one can determine each topological 
amplitude from experimental data with the help of fitting, if a sufficiently large number 
of measurements are available. So far several works have been done to study hadronic 
charmed meson decays in the framework of the quark diagram approach. However, because 
of the difficulty to manage the large number of parameters in fitting, those works have used 
at least in part certain model-dependent information {e.g., information on SU(3) breaking 
from factorization or QCDF, etc.) to obtain a £t. 

In this work we shall study direct CP asymmetries of SCS D(s) —)■ PP (P = tt, iP, 
decays using the updated experimental data. For a least model-dependent analysis of the 
charmed meson decays, we choose the quark diagram approach and perform the analysis 
in the most general way, i.e., without using any model-dependent information The present 
experimental data show that the measured values of branching fractions of —)■ PP 

are quite accurate, but direct CP asymmetries include very large errors. In this situation 
our strategy is: (i) first, to fit the experimental data, especially the accurately measured 
branching fractions, by using the topological amplitudes as parameters, (ii) then, to extract 


^ This will be Case I of our analysis. In Case II we shall impose certain constraints on two of the parameters 
encouraged by the fit obtained for CF charmed meson decays. 
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each topological amplitude from the ht, and (iii) to make predictions for the direct CP 
asymmetries. Since the large number of parameters are involved in this analysis, obtaining 
a satisfactory £t without using any mo del-dependent information turns out to be a very 
difficult task. 

We would like to emphasize that (as we shall see later) for direct CP asymmetries of SCS 
D decays, the penguin amplitudes corresponding to the internal 6-quark [i.e., internal b- 
quark penguin amplitudes) multiplied by the CKM factor V^^yub become important so that it 
cannot be ignored. In our analysis we shall explicitly express all the relevant 6-quark penguin 
amplitudes, including the 6-penguin exchange and the 6-penguin annihilation ones. Then, 
we shall determine the magnitudes and strong phases of these 6-quark penguin amplitudes. 
This is one of different points of our work from the previous other works. In other words, all 
possible strong penguin contributions including the internal 6-quark penguins are explicitly 
included in our analysis. We shall see that the internal 6-quark penguin contributions are 
comparable with the internal s- and d-quark penguin ones. 

This paper is organized as follows. In Sec. II, we introduce the Wolfenstein parametriza- 
tion of the CKM matrix up to order A® and topological quark diagrams relevant to iA —)■ PP 
decays. In Sec. Ill, the explicit SU(3) decomposition of the decay amplitudes and its rele¬ 
vance to direct CP asymmetries are presented. In Sec. IV, we perform the analysis by 
taking into account SU(3) breaking to determine the topological amplitudes, and predict 
the direct CP asymmetries. Our conclusions are given in Sec. V. 


II. FRAMEWORK 


It is well known that the 
obtained from the standard 
up to order A® are given by 


Wolfenstein parametrization of the 
Chau-Keung(CK) parametrization. 


43| 


CKM matrix can be easily 
The CKM matrix elements 
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2 8 

~T 6 ’ 


A , 


AX^{p — if]) 


\ 


V, 


(CK) 


Wolf 


-X+^A‘^{l-2p-2i-n) , 1-^-^(1 + 4^2) 

- ^[l-4^^(l-4p-4i77)] , 




+ 0(A") . 


AlA^(l — p — ip) 
-\-^A{p + ip) , 


-^A2 + ^A(l -2p- 2ip) 1 - ^^2 

+ ^Al, 


( 1 ) 


Regarding charm decays, the relevant CKM factors are Xg = V*^Vuq {q = d, s, b). From the 
above matrix elements, one hnds 

Xd = —A + — + —(1 + AA^) — X^A‘^{p + irj) + 0{X'^) = X^J^ + A® + O(A^), 

As = A— - -^(1 + 4A^) + C>(A^), 

A, = X^A\p - ip) + 0{A) = |A,|e-*^ + 0{A), (2) 


and A® = —X^A^{p + ip) = = “|Ab|e+*'^ 


= _| \(2)|p+i7 = 


where A^^^ = -A + f + A^ [l + Xa^ 

Notice that the imaginary terms appear in both Xd and X;, at order A® Thus, for CP 
asymmetries in charm decays, the internal 6-quark contributions (penguin contributions) 
with the CKM factor X/, as well as the d-quark ones with A^ become important. 

In SU(3) decomposition of the decay amplitudes for D —)■ PP [P = vr, iC, p^'^) modes, the 
decay amplitudes are expressed in terms of topological quark diagram contributions. The 
topological amplitudes corresponding to different topological quark diagrams, as shown in 
Figure 1, can be classihed into three distinct groups as follows: 


1. tree and penguin amplitudes: 

• T, color-allowed tree amplitude (equivalently, external hF-emission); 

• C, color-suppressed tree amplitude (equivalently, internal hF-emission); 

• V, QCD-penguin amplitude; 


2 Note that up to order of A®, Xd —X — A®e+®'>', Ag ^ A, and Xb ^ X^e Therefore, ~ A^ ~ 

I As(d) \ 

2 X fO-h 
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• S, singlet QCD-penguin amplitude involving SU(3)-singlet mesons (e.g., 

r/W, uj, (/)); 

• 'Pew ■ color-favored EW-penguin amplitude; 

• Pew ■ color-suppressed EW-penguin amplitude; 

2. weak annihilation amplitudes: 


• S, lE-exchange amplitude; 

• A, lE-annihilation amplitude; {E and A are often jointly called “weak annihila¬ 
tion” amplitudes.) 

• P£, QCD-penguin exchange amplitude; 

• PA, QCD-penguin annihilation amplitude; 

• P£ew '■ EW-penguin exchange amplitude; 

• P^EW : EW-penguin annihilation amplitude; {PE and PA are also jointly called 
“weak penguin annihilation”.) 

3. flavor-singlet weak annihilation amplitudes: all involving SU(3)F-singlet mesons, 


• SS, singlet lE-exchange amplitude; 

• singlet lE-annihilation amplitude; 

• SP8, singlet QCD-penguin exchange amplitude; 

• SPA, singlet QCD-penguin annihilation amplitude; 


SPE-ewi '■ singlet EW-penguin exchange amplitude; 
SPAew '■ singlet EW-penguin annihilation amplitude. 


The reader is referred to Ref. 


44l | for details. Each topological quark diagram in this 


approach includes all possible strong interactions to all orders. Therefore, analyses of topo¬ 
logical diagrams can provide information on final-state interactions. 


III. DECAY AMPLITUDES AND CP ASYMMETRIES 

The decay amplitudes of two-body hadronic D decays can be represented in terms of the 
topological quark diagram contributions. In general the decay amplitudes of 
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FIG. 1: Topology of possible flavor diagrams (not Feynman diagrams: Each topological diagram 
includes all possible strong interactions to all orders.): (a) color-allowed tree T, (b) color-suppressed 
tree C, (c) QCD-penguin V, (d) singlet QCD-penguin S diagrams, (e) VF-exchange £, (f) W- 
annihilation A, (g) QCD-penguin exchange VS, and (h) QCD-penguin annihilation VA diagrams. 
The color-suppressed EW-penguin color-favored EW-penguin PeW; EW-penguin exchange 

VSew and EW-penguin annihilation VAew diagrams are obtained from proper replacement of 
gluon lines by Z-boson or photon lines in (c), (d), (g), (h), respectively. 


and ^ K+R- can be written as 


A{D° ^ TT+TT-) = Xd{r + S + Vd + VAd + V£dY^ + XsiVs + VAs + VEsT^ 

+ xYv, + vA, + vE,y^- 

A{D^ ^ K^K-) = xYr + E + Vs + vAs + + xyVd + VAd + VEdf^ 


Xbi'Pb + T^-^b + T^Ei 


\KK 


( 3 ) 


( 4 ) 


where Vg, VAg, VEg are the QCD-penguin, QCD-penguin annihilation, QCD-penguin ex¬ 
change amplitudes with an internal g-quark, respectively. Here the internal 6-quark contri¬ 
butions with the CKM factor Xb are explicitly shown. ^ 

The possible sources of CP violation in charm transitions include: CP violation in AC* = 1 
decay amplitudes (direct CP violation) and CP violation through AC = 2 mixing 

(Indirect CP violation). The latter has been estimated to be zero in recent experiments 241. 


^ For simplicity, electroweak-penguin contributions have been neglected since they are expected to be very 
small, as previous works in the literature. 
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The direct CP asymmetry is defined as 


o>cp{D —^ /) — 


B{D ^ /) - B{D ^ /) 


(5) 


B{D^f) + B{D^fy 
where B denotes a branching fraction. For D ^ f decay having two strong and two weak 
phases, the decay amplitude can be generally written as 


A{D f) — \Ai\e 


-iSi -i<pi 


+ I A 2 1 e 


—*<52p—*02 

C' 1 


( 6 ) 


where 61^2 and are strong and weak phases, respectively. Then, the direct CP asymmetry 
is given by 


o>cp{D —^ /) — 


2 Ml AI 2 

sin(( 5 i - 

- 62) sin((^i - (j)2) 

lAI 

P + l^2| 

p -|- 2 Ml M 2 

COs( 5 i - 62) COs((/)i - (j)2) ’ 


(7) 


which is non-zero only if A(5 = (5i — (52 7 ^ 0 , A(/) = (/>i — (/>2 7 ^ 0 , and |Ai^ 2 | 7 ^ 0 . 

To study the direct CP asymmetry, the decay amplitude of —)■ K^K~, for instance, 

can rewritten as 


AiD^ ^ TT+TT-) = Ar + A^, 


( 8 ) 


where up to order A®, 

= Xy\T + + PEdP^^^) 

+ \s{PsP^^ + PAsP^^^ + PEsP^^^) 

= \Ar\e-^^\ 

= xf\T + EP^^ + PdP^^ + PAdP^^^ + PEdP^^^) 

+XpPi,P^^ + PAbP^p^ + PEbP^p^) 

= -|Ai^^|(T + EP^^ + PdP^^ + PAdP^^^ + PEdP^^^)e+P 
+ \Xb\{PbP^^ + PAbP^p^ + PEbP^^^)e-^'^ 

= |^2l|e"*'^'e+*^ +(9) 

explicitly with the strong phases (5’s and the weak phase 7 , such as = Ee^^^. Here 
T(= Te*^^) is taken to be real (he., 5t = 0) and the other strong phases (5’s are the relative 
ones to 6 t- We note that the weak phase 7 appears only in the terms having the CKM factor 
Aft or A^^^ of order A^. Thus, one expects that the direct CP asymmetry for —)■ 7r'''7r is 
CKM-suppressed by the factor of A^ .V. 10-3. 

Within the SM, for CP asymmetries in charm decays one has to go to the CKM matrix 
through order A® in the Wolfenstein parametrization 


45l | and understand the differences 





TABLE I: Topological amplitudes for D —?■ PP {P = vr, K, modes (singly Cabbibbo-suppressed 
decays). The strong phases 5’s and the weak phase 7 are explicitly denoted. The CKM factor 
Ad = where A^^^ = —A + ^ + ^(1 + 4A^) and A^^^ = —\^A^{p + ir]) = — |A^^^|= 

— I Able"*'*'’'. For simplicity, singlet QCD penguin, flavor-singlet weak annihilation and electroweak- 
penguin amplitudes have been neglected. 


Mode 

Representation 

D —>• TT^TT” 

- |A^^^|e+*T')(r -k -k Pde^^^ + PA^P^p^ + PEdP^P^) 

+Xs{PsP^^ + PAsP^P^ + PEsP^Pp) 

0 

0 

-k|Ab|(Pfee*'^^’ -k PAbP^P'^ + PEbP^pp)e-P 
-^[{X^d^ - \xf\e+P){-CP^c + EP^p + PdP^p + PAdP^p^ + PEdP^pp) 
+Xs{PsP^p + PAsP^P^ + PEsP^pp) 

TT^Vr® 

+|A6|(Pbe*^^ + PAbP^P^ + PEbP^p^)e-P] 

K+K- 

i^d^ - \ ^d^\eP'^){PdP^P + PAdP^P^ + PEdP^PP) 

-kAs(r -k EP^p + PsP^P + PAsP^P^ + PEsP^Pp) 


-klAfeKPfoe^'^^ -k PAbP^p^ + PEbP^pp)e-P 
i^d^ - \X^d^\e+P){EP^P + 2PAdP^P^) 

+Xs{EP^p +2PAsP^P^) 

K+K^ 

+\Xb\{2PAbP^p^)e-P 

i^d^ - \X^d^\e+P){AP^^ + PdP^P + PEdP^pp) 

+Xs{T + PsP^P + PEsP^PP) 

+ \Xb\{PbP^P + PEbP^pp)e-^'^ 

TT^T] 

i^d^ - \^d^\e-^^^){-EP^^ + PdP^P + PEdP^PP) cos (1) 

—Xs[-^C sin (l)p^<p -k {PsP'^P + PEsP^Pp) cos 0)] 

+1 AftI ) cos (/>e-*T' 

TT^T]' 

7r~^r] 

- \XP\e+E)[-EP^P + PdP^P + PEdP^PP) sin0 
i^P - + CP^c + 2AP^^) + V2{PdP^P + PEdP^PP)] cos 0 

-kAs[—C sin (jyp-^^^ + ^/2{PsP^P + PEgP^Pp) cos 4>] 

TT~^ri' 

+\Xb\[V2{PbP^P + PEbP^PP)] cos (/>e-*T' 

(^P - + CP^c + 2Ae*^^) + V2{PdP^P + PEdP^PP)] sin^ 

+ As[C' cos 4>P^^ -k \/2{PsP'^P + PEsP^Pp)s\n(l)\ 

+ \Xb\[V^{PbP^P + PEbP^Pp)] sm.<j)e~^'^ 
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TABLE II; {Continued from Table I) Topological amplitudes for Dg —)• PP {P = TT,K,r]^'^) modes 
(singly Cabbibbo-suppressed decays). 


Mode Representation 

Dg tt+ATO - |A^^V+*T')(r + Pae^^p + PEae^^pp) 

+Xg{Ae^^^ + PgC^^P + PEge^^Pp) 

+ |A6|(Pfee*^^ + PEbe^^PP)e-^'i 

K+tt^ ^ - |Aj^|e+*T')(-Ce*^c + p^e^Sp ^ pe^P^pe) 

+Xg{Ae^^^ + Pge^^P + PEge^^Pp) 

+\Xb\{Pbe^^P + PEbe^^pp)e-^^ 

K+r] (A^^^ - Ile+^T-)+ p^pi^P + pp^e^^pp ) cos </> - {Pae^^P + PE^e^^pp ) sin 4>] 

+As[;^(Ae*'^"^ + PdP'^P + PEdC^^Pp) coscj) 

-{T + C'e^'^c- + Ae^^^ + PgC^^P + PEgC^^Pp ) sin (p] 
+|Afe|(Pfee*'^^ + PEbe^^Pp){^ coscf) - sm(f))e~^'^ 

K+rj' (A^^^ - |Ai^^|e+*T')[^(Ce*^c + Pde^^P + ) sin <(. + {Pde^^P + PEde^^PP)cos 4>] 

Xs[^{Ae'‘^^ + PdP^^ + ) sin (f 

+(T + Ce^^'P + + PgC^^P + PEgC^^Pp) cos (j)] 

lAfeKPfce^-^^ +PPfee^-^^^)(^ sincj) + cos(p)e 


between other parameterizations 43|. The SM generates CP asymmetries of order A^ for 
SCS charm decays and does not generate any CP ones for DCS decays np to order of A®. 

In Table I and II, the decay amplitndes of D —)■ PP (P = vr, JP, modes are expressed 
in terms of the topological amplitudes explicitly with the strong phases (5’s and the weak 
phase 7 . For D —)■ and Dg Krj^'^ modes, the rj — rj' mixing is considered: in the 

flavor basis. 


COS( 

sin ( 


— sm ( 
cos 6 


Vq 

.Vs. 


( 10 ) 


= 40.4° 


46|. 


where Vg = ^ + dd ), rig = ss, and the mixing angle 0 = 

In order to understand dynamics behind CP violation, one hrst needs to obtain decay 
amplitudes of the corresponding decays. The only model-independent way to analyze these 
amplitudes is to consider contributions from all possible quark diagram processes. We would 
like to point out that unlike all the previous works, we consider contributions from all pos¬ 
sible strong penguin diagrams including 6 -quark processes. However, limited experimental 
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TABLE 

decays). 


Ill: Experimental data for D —)■ PP {P = 7r,K,ri^'^) modes (singly Cabbibbo-suppressed 
Branching fractions and direct CP asymmetries are shown in units of 10“^ D- 


Mode Expt Br (xlO 


—)• 7r'*‘7r“ 

—)■ TT^TT^ 

—> TT^TT^ 

K+K- 
D+ K^K+ 

—)• TT^T] 

—>■ TT^T]' 

—>■ TT^T] 

D+ -5- TT+T/' 
D+ K^7r+ 
D+ K+tt^ 
D+ K+t] 
D+ K+i 


1.400 

± 

0.026 

1.19 

± 

0.05 

0.82 

± 

0.05 

3.96 

± 

0.08 

5.66 

± 

0.32 

0.34 

± 

0.08 

0.68 

± 

0.07 

0.89 

± 

0.14 

3.53 

± 

0.21 

4.67 

± 

0.29 

2.42 

± 

0.16 

0.63 

± 

0.21 

1.75 

± 

0.35 

1.80 

± 

0.60 


Expt acp (xlO 
2.2 ± 2.1 
29.0 ±29.0 
0.0 ±50.0 
-2.1 ± 1.7 
-1.1 ±2.5 
230 ± 190 


10.0 ± 15.0 
-5.0 ± 12.0 
12.0 ± 10.0 
-266 ± 238 
93.0 ± 152 
60.0 ± 189.0 


information and the large nnmber of parameters to be determined has made it a rather 
difficult task. 

IV. THE analysis AND SU(3) BREAKING EFFECTS 

The topological amplitudes can be extracted from the available experimental information, 
such as branching fractions and direct CP asymmetries of SCS charm decays. Previous 
studies show that it is difficult to £t the experimental data in the x^ analysis without using 
additional theoretical information such as QCDF and pQCD. It is because the number 
of parameters in theory are large and the experimental information particularly on acp 
measurements in charm sector is rather poor. Also, the SU(3) breaking effects cannot be 
ignored, as it is impossible to ht the data within flavor SU(3) symmetry. 
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A. The analysis in SU(3) limit 


We perform the analysis with 25 observables (experimental branching fractions and 
direct CP asymmetries) as inputs, shown in Table HTTl Within SU(3) symmetry, the total 
number of parameters is 20 , including the magnitudes and strong phases of topological 
amplitudes, as shown in Table HVl Although the degree of freedom (d.o.f.) in this case is 5, 
it turns out that it is impossible to obtain a £t with an acceptable x^/d-o.f. In this fit, we 
hnd 650 < y^/d.o./. < 1000, being unacceptably large. This fact has also been supported 
by previous works based on a similar analysis [7|. 

B. The analysis including SU(3) breaking effects 

In order to include SU(3) breaking effects, we divide decay processes in categories of 
TTTT, KK, Ktt and Krj^''' for D and Dg decays. Then we introduce hve additional 

parameters in their amplitudes, namely, ^k-k and where each App 

characterizes each category of PP modes. The SU(3) broken decay amplitude in such case 
may, for example for —)■ tt+tt", be given as 

^ TT+TT-) = - |Ai^^|e+*^)(T + + P^P^^ + PA^P^^^ + PE^P^^^) 

+ XsiPgP^^ + PAgP^p^ + PEgP^p^) 

+ \M{PbP'^ + PA,P^^^ + PEhP'^^)e-^^ ](1 + A,^). (11) 

It may be noted that introduction of SU(3) breaking parameters lowers the x^/d.o./. to 
acceptable limits. To perform the x^ analysis, we use 25 observables (experimental branching 
fractions and direct CP asymmetries) as inputs, shown in Table IIIII The total maximum 
number of parameters in this analysis is 24, including the magnitudes and strong phases of 
topological amplitudes and SU(3) breaking parameters, as shown in Table HVl Based on the 
choice of parameters we have two cases as follows. 

1. Case I: fit without any constraints 

We perform a complete analysis without using any constraints to extract 24 parameters 
with the weak phase 7 = 63°. In this case the degree of freedom is one. The obtained 

^ The CP asymmetry of —>■ has not been used as an input in the fit. 
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TABLE IV: Case I: Topological amplitudes and SU(3) breaking parameters determined from the 
fit without any constraints. The magnitudes and strong phases of the amplitudes are given in units 
of 10“® GeV and degrees, respectively. 


No. 

Parameter 

Value 

No. 

Parameter 

Value 

No. 

Parameter 

Value 

1 

T 

3.69 

9 

PAg 

1.29 

17 


-144.10 

2 

C 

2.91 

10 

PAb 

1.47 

18 

SpA 

92.22 

3 

A 

1.50 

11 

PEd 

0.30 

19 

SpE 

24.50 

4 

E 

1.65 

12 

PEg 

2.46 

20 

ZXtj-TP 

0.006 

5 

Pd 

2.18 

13 

PEb 

0.60 

21 

Ai-K-K 

0.07 

6 

Ps 

1.00 

14 

Sc 

145.15 

22 

57 

<! 

0.21 

7 

Pb 

1.35 

15 

Sa 

-143.53 

23 

^KK 

0.011 

8 

PAd 

0.32 

16 

Se 

-110.73 

24 


0.17 


numerical values of all the parameters with x^/d-O-f- = 8.0 are listed in Table HV] For 
instance, it is found that 


T = 3.69, 


Pd = 


( 12 ) 


Interestingly, the htted values for T, C, E are quite close to the values obtained by other 


works [ 7 , for analysis of Cabibbo-favored modes of charm decays, except for the sign of 

the phases It is noted that topological amplitudes for the internal 6-quark penguin are of 
the same order as the s- and d-quark penguin amplitudes, though Pg ^ Ph < Pd 


® They^/d.o./. value does not seem small. However, this feature has been known in the previous works 
101123 using the similar global fit to H —?► PP data. The y^/d.o./. value in the present fit mainly comes 


from the differences between theory fits and the central values of the poorly measured data for direct 
CP asymmetries. For the CP asymmetries, several modes in this fit induce sizable values: e.g., for 
—>■ 7r+7r“, 7r°7r+, K~^K~, Dg —>■ A'°7r+, A+7r°, the corresponding values are 0.8, 1.0, 0.87, 

1.4, 1.2, 0.8, respectively. 

® In our analysis we have carefully examined the value of y^/d.o./. to determine the magnitudes and phases 
of the amplitudes. If the values of the phases change from the above ones, the y^/d.o./. increases. 

^ It may be argued that contributions from penguin amplitudes are usually expected to be considerably 
smaller in magnitude as compared to the color-allowed tree contribution. However, they may be enhanced 
due to long-distance FSI so as not to be neglected. 
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For D —)■ 7171, KK and Dg —?■ Ktt modes, the SU(3) breaking effects are found to be only 
a few percents or less. For D —)■ and Dg modes, the SU(3) breaking becomes 

as large as about 20 %. 

Using the htted parameter values, the branching fractions and direct CP asymmetries 
are predicted, as shown in Table |Vl We find that the £t obtained for the branching fractions 
is surprisingly good (Compare Table IVl with Table HTTl ). Almost all the predicted branching 
fractions are in very good agreement with the experimental data. It should be emphasized 
that the branching fractions of tt+tt" and K^K~, having been a long standing 

puzzle, are also in excellent agreement with the data. The only ambiguities in fitting arise 
for the a^p’s of the observed decay modes which seems not to fit well in the present scenario. 
Since the experimental values of CP asymmetries still have very large uncertainties, in order 
to make a reliable conclusion, one has to wait for more precise CP measurements from the 
future experiments such as LHCb and Belle II. We would like to remark that to the lowest 
order acp{D PP) will include the 6 -quark penguin contribution multiplied by sin 7 (which 
is close to 1 for 7 = 63°). Thus, the direct CP asymmetries of D PP are expected to get 
a signihcant effect from the sizable contribution of the 6 -quark penguin shown in Table IIVI 
From the fitted parameters, one can make pure predictions for direct CP asymmetries of 
—)■ K^K^, 7 r °7 and modes which have not been used as the inputs. They are found 
to be 


acp(T>° ^ K^K^) = 1.42 x 10-^ 
acp{D^ —)■ 71^7]) = 2.8 X 10“'^, 

acp(T>° ^ A') = -6.8 X 10"^ (13) 

The direct CP asymmetry for —)■ K^t] is expected to be large, being of 0(10“^). For 

—)■ and —)■ tt^t] modes, a^p’s are expected to be of (P(10“^). For —)■ 

TT^r]', D~^ —)■ and —)■ modes, acp’s are close to 0(10“^). For the other modes 

including —)> tt+tt", acp's are of 0(10“^) or smaller. Future experimental measurements 
for these observables will provide useful information on this scenario. 

2. Case II: fit with C = 0.8T and Pg = (1.0 ± 0.2)Pfi 

We put some constraints on our analysis encouraged by the fit obtained by for Cabibbo- 
favored decay modes of charmed mesons. With the constraints C = 0.8T and Pg = (1.0 ± 
0.2)Prf, we obtain the fit for the remaining 22 parameters. The relation C = 0.8T holds 
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TABLE V; Case I : Branching fractions and direct CP asymmetries (in units of 10 obtained 
from the ht without any constraints. 


Mode 

Br (xlO 

o-cp (xlO 

—>• 7r'’'7r“ 

1.40 

0.30 

—>■ 7r^7r+ 

1.19 

0 

TT^TT^ 

0.83 

0.86 

K+K- 

3.96 

-0.51 

D+ K^K+ 

5.68 

-0.76 


0.36 

1.42 

> TT^ry 

0.68 

0.28 

—>■ vr'^Ty' 

0.89 

-0.68 

—)■ T^^rj 

3.50 

1.50 

—>• 7r'’“r/' 

4.65 

-1.02 

D+ K^-k+ 

2.41 

0.30 

D+ A+ttO 

0.66 

-0.93 

Dt K+r] 

1.56 

-42.0 

D+ K+r]' 

2.22 

0.47 


TABLE VI: Case II: Topological amplitudes and SU(3) breaking parameters determined from the 
fit with C = 0.8T and Pg = (1.0 ± 0.2)Pd. The magnitudes and strong phases of the amplitudes 
are given in units of 10“® GeV and degrees, respectively. 


No. 

Parameter 

Value 

No. 

Parameter 

Value 

No. 

Parameter 

Value 

1 

T 

3.13 

9 

PAg 

1.01 

17 

5p 

-145.11 

2 

C 

2.50 

10 

PAb 

1.48 

18 

5pA 

81.46 

3 

A 

1.26 

11 

PEd 

0.52 

19 

5pB 

21.65 

4 

E 

1.96 

12 

PEg 

2.42 

20 

^TTTT 

0.003 

5 

Pd 

1.84 

13 

PEb 

0.50 

21 

^Ktv 

0.35 

6 

Ps 

1.84 

14 

5c 

137.54 

22 

<! 

0.35 

7 

Pb 

1.50 

15 

5a 

-159.66 

23 

Abb 

0.33 

8 

PAd 

0.92 

16 

5b 

-113.50 

24 

^Kt] 

0.32 
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also in the result of Case I. In this case the degree of freedom is 3. The htted values of the 
parameters with x^/d-o.f. = 9.6 are shown in Table IVTl Though the parameter values are 
similar to those obtained in Case I, in this case Pb ^ Pd = Ps, and the SU(3) breaking effects 
are as large as 32 ~ 35% for all the modes, except tttt modes where the SU(3) breaking turns 
out to be very small. 

Table IVIII shows the predicted branching fractions and ocp’s. Here also the fit for 
branching fractions is excellent. As in Case I, the branching fractions of tt+tt" 

and K+R- are also in excellent agreement with the experimental data. In this case, 

compared with Case I, more pure predictions for direct CP asymmetries can be made. From 
the htted values of the parameters, we hnd for —)■ tt+tt", K~^K~, K^K^, ir^rj and Tr^rj' 
modes which have not been used as the inputs, 

acp{D^ = 3.3 x 10““^, 

acp(T>° ^ K+R-) = -7.3 x 10-^ 
acp(T>° ^ R^R^) = -5.8 x 10-^ 
acp{D^ TT^r]) = —4.7 x 10“'^, 

acp{D^ TT^V') = -3.3 X 10"^ (14) 

In this case, the direct CP asymmetries for several modes are expected to be smaller, com¬ 
pared with Case I. The direct CP asymmetries for —)■ R'^n^ and R^rj are of (P(10“^). 

For —)■ TT^Ti^, R~^R~, ®cp Eire close to 0(10“^). For the other modes 

including tt+tt", acp’s are expected to be of 0(10“'^) or smaller. Future experiments 

such as upgraded LHCb and Belle II will help to determine which scenario is more reliable. 


V. CONCLUSIONS 

We have performed a model-independent analysis of CP violation in the singly Cabibbo- 
suppressed D —)■ PP decays. In light of the most recent CP asymmetry measurements a^p 
for the —)■ R^R~ and —)■ 7r'''7r“ modes, we have analyzed direct CP violation in 
D —)■ PP decays within the SM. It is believed that direct CP violation in such decays may 
arise from the interference between various topological amplitudes and the asymmetries are 
expected to be of (P(10“^) or smaller. We have taken advantage of the quark diagram ap¬ 
proach to extract the topological amplitudes and relative strong phases from the measured 
experimental data. We have also taken care of SU(3) breaking effects which are considered 
to be large in these decays. It may be noted that for charmed meson decays it is difficult 
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TABLE VII: Case II : Branching fractions and direct CP asymmetries (in units of 10 obtained 
from the fit with C = 0.8T and Pg = (1.0 ± 0.2)Pd. 


Mode 

Br (xlO 

acp (xl0"3) 

—>• 7r'’'7r“ 

1.40 

0.33 

7r®7r+ 

1.18 

0 

—>• 7r®7r® 

0.89 

0.71 

^ K+K- 

3.98 

-0.73 

D+ K^K+ 

5.56 

-0.17 

D® ^ A® A® 

0.39 

-0.58 

—>• TT^T] 

0.68 

-0.47 

—> TT^g' 

0.91 

-0.33 

—>■ 7r~^g 

3.39 

0.89 

> TT 

4.57 

-0.67 

D+ 

2.41 

-0.15 

D+ A+vr® 

0.54 

-2.43 

D+ K+g 

1.48 

1.4 

D+ K+g' 

2.32 

-0.23 


to induce CP violation at tree level. Thus, one has to include contributions from penguin 
amplitudes at certain level. In order to explain the experimental branching fractions, we 
have included contributions from all possible strong penguin amplitudes involving penguin 
exchange and penguin annihilation diagrams. We have already stated that due to the large 
uncertainties in measured acp’s, CP violation in charmed meson decays has not yet estab¬ 
lished experimentally. This indicates that a more precise approach is needed to analyze 
CP violation in these decays. Therefore, we have included all the possible internal 6-quark 
contributions to penguin diagrams (up to order of A®) which have been ignored so far in the 
previous works. 

In the quark diagram approach, we have determined topological amplitudes for D —?■ PP 
decays through the analysis by using the available experimental branching fractions and 
direct CP asymmetries. Our analysis has been performed in the most general way, i.e., 
without using any model-dependent informationon the parameters (Case I). We have divided 
our analysis into two cases: without any constraints (Case I) and with certain constraints 
(Case II), having the best y^/d.o./. = 8.0 and 9.6, respectively. Consequently, we have 
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predicted the unmeasured direct CP asymmetries. Our findings are summarized as follows. 

1. We observe an excellent fit for branching fractions, including B{D^ —)■ tt+tt”) and 

—)■ K~^K~), in both the cases as shown in tables IVl and IVIIl 

2. It may be argued that contributions from penguin amplitudes are considerably smaller 
in magnitude as compared to the color-allowed tree contribution. However, they may 
be enhanced due to long-distance FSI resonances so as not to be neglected. 

3. We find that penguin contributions from the internal &-quark can be sizable and es¬ 
pecially for CP asymmetries in charm decays, they become important and cannot be 
neglected. 

4. Inclusion of 6-quark contributions result in CP asymmetries that lies in range within 
the SM: i.e., 10“'^ < acp < 10“^ for both the cases. 

5. In Case I: (i) The observed T, C and S acquire values closer to those found in a similar 
analysis based on Cabibbo-favored charm decays. The SU(3) breaking is up to 20% 
for TTrj^''^ and Kr]^'^ modes. 

(ii) We predict acps for /ttt]/ttt]' modes to be of 0(10“^) ~ 0(10“'^). 

6. In Case II: (i) We use constraints C = 0.8T and Pg = (1.0 ± 0.2)Prf to increase 
our predictability and to test our fit. We predict acp{T^^T^~) = 0.33 x 10“^ and 
acp{K^K~) = —0.73 x 10“^ with Aacp = —1.00 x 10“^, which are consistent with 
the recent experimental result. 

(ii) We also predict acp^s for /7rr]/7rr]' modes to be of the same order of magnitude 

and sign. 

(iii) SU(3) symmetry breaking as large as 35% is observed in Case II for the modes 
other than tttt. 


We wish to remark here that our results are more or less comparable/consistent with other 
theoretical works based on similar approaches. In order to get a clearer picture, more 
precise experimental measurements are required to determine the exact magnitude of CP 
asymmetries in charmed meson decays. We hope that such an analysis would be helpful in 
diagnosing possible evidence of new physics in the charm sector. 


* Note added: When finishing the paper, we have found that two new works 22, l23(] have 
just come out. Although those two works also use the similar approach to ours, there are 
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clear differences between our work and theirs. As mentioned in Introduction and Conclu¬ 
sions, we have performed the analysis in the most general way, i.e., without using any 
model-dependent information on the parameters. In order to get more precise and reliable 
estimates of the amplitudes, we have included contributions from all possible strong penguin 
amplitudes, including the internal 6-quark penguins. Then, we have made predictions for 
direct CP asymmetries of SCS D decay modes, including yet unmeasured ones. 
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